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Abstract 

We study the deformation kinetics for the (glassy) sohd hehum confined in microscopic domain 
at very low temperature regime by using an absolute-reaction-rate model considering the shear 
thinning behavior which means, once material being subjected to high shear rates, the viscosity 
diminishes with increasing shear rate. Our calculations show that there might be nearly friction- 
less fields for rate of deformation due to the almost vanishing shear stress in microtubes at very 
low temperature regime together with wavy-rough corrugations along micropores and the slip. 
As the tube size decreases, the surface-to-volume ratio increases and therefore, surface roughness 
will greatly affect the deformation kinetics in micropores. After using the boundary perturbation 
method, we have obtained a class of temperature and activation energy dependent fields for the de- 
formation kinetics at low temperature regime with the presumed small wavy roughness distributed 
along the wall of an cylindrical microtube. The critical deformation kinetics of the glassy matter is 
dependent upon the temperature, activation energy, activation volume, orientation dependent and 
is proportional to the (referenced) shear rate, the slip length, the amplitude and the orientation of 
the wavy-roughness. Finally, we also discuss the quantitative similarity between our results with 
Ray and Hallock [Phys. Rev. Lett. 100, 235301 (2008)]. 

PACS numbers: 67.80.K-, 67.80.B-, 67.25.dr, 83. 60. St, 83.60.Rs,83.50.Lh 
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I. INTRODUCTION 



Grigor'ev et al. recently performed high-precision pressure measurement in sohd ^He 
samples grown by capillary blocking technique. In all their nonannealed hep crystals, the 
temperature dependence of pressure demonstrates a contribution proportional to T^, the 
latter becomes the leading term at temperatures T < 300 mK, at which supersolid effects 
were observed. Grigor'ev et al. thus claimed that such a behavior may be ascribed to a 
glassy phase. They also found that this glassy contribution to pressure can be eliminated 
only by a substantial annealing : A dramatic pressure decrease of ~ 2 bar was observed 
under annealing at temperatures very close to the melting point. They thus conjectured 
that this effect is due to sohdification of liquid or glass captured in closed cavities during 
the growth process [1]. 

Meanwhile Day and Beamish recently observed an approximately 10% increase in the shear 
modulus of "^He at low temperatures (below 200 mK) [2] . Ray and Hallock found evidence 
for flow through solid ^He in at least some cases [3]. To be specific, Ray and Hallock have 
conducted experiments that show the first evidence for flow of helium through a region 
containing solid hep ^He off the melting curve. Their phase diagram appears to have two 
regions. Samples grown at lower pressures show flow, with flow apparently dependent on 
sample history, with reduced flow for samples at higher temperature, which is evidence for 
dependence on temperature. Samples grown at higher pressures show no clear evidence for 
any such flow for times longer than 10 hours [3]. The temperatures utilized for their work 
are well above the temperatures at which much attention has been focused, but interesting 
behavior was seen. 

Note also that Clark et al. have studied the thermal history of the resonant frequency of a 
torsional oscillator containing solid ^Hc. They found that the magnitude of the frequency 
shift that occurs below ~ 100 mK is multi-valued in the low-temperature limit, depending 
strongly on how the state is prepared [4]. However, Anderson interprets the observed NCRI 
(nonclassical rotational inertia, please see [4] for the detailed references) as a consequence 
of vortex liquid and supercurrents flowing around a thermally excited fluctuation of vortices 
[5]. The subsequent experiments [6] reported the flrst observations on the time- dependent 
dissipation when the drive level of a torsional oscillator containing solid "'He was abruptly 
changed. The relaxation of dissipation in solid ''He shows rich dynamical behavior including 
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exponential and logarithmic time-dependent decays, hysteresis, and memory effects. As re- 
marked in [6] , their procedure of initiating the oscillation at a high temperature is likely to 
have brought the sample solid ^He to a quasi-equilibrium array of vortices for the drive level 
at the temperature of measurement. The observed logarithmic time dependence probably 
arises from the therm,ally activated motion of vortices from one metastable state to another. 
This process continues throughout the measurement while the sample was maintained at the 
constant temperature. Just after the drive level has been decreased, the existing vortices are 
out of equilibrium with the new drive level. The observed undershoot may be a transient 
excess dissipation owing to the 'wrong' number of vortices. The vortices then must adjust 
to the new drive level by, say, moving out of the sample. The characteristic time involved 
in the macroscopic motion of the vortices was then the observed time constant [6]. The 
motion presumably involves processes occurring both within the sample and at the surface 
boundaries. The time constant would become shorter as the vortices begin to move more 
freely. The proposed vortex liquid state [5] appears to occur above about 60 mK. When the 
vortex liquid 'freezes', the time constant would diverge [6]. 

Researches mentioned above all imply that the supersolid (material) is not so simple but 
really complex. We also noticed that as reported by Dyre et al. in 1996 for glass-forming 
molecular liquids the (high-frequency) shear modulus increases as the temperature decreases 
[7]. The latter resembles that reported in [2] by Day and Beamish. Thus, it's necessary to 
study the deformation kinetics as well as transport of amorphous and/or glassy material 
(presumed solid helium to be almost the same) under confined microdomains at rather low 
temperature regime! 

Glasses arc amorphous materials of polymeric, metallic, inorganic or organic type. The 
plastic deformation of amorphous materials and glasses at low temperatures and high strain 
rates is known to be inhomogeneous and rate-dependent. In fact, the mechanical behavior 
of amorphous materials such as bulk metaUic glasses [7-10] continues to present great theo- 
retical challenges. While dislocations have long been recognized as playing a central role in 
plasticity of crystalline systems, no counterpart is easily identifiable in disordered matter. 
In addition, yield and deformation kinetics [11-12] occur very far from equilibrium, where 
the state of the system may have a complex history dependence. 

In recent years, considerable effort was geared towards understanding how glasses respond 
to shear [12]. Phenomena such as shear thinning and 'rejuvenation' are common when shear 
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deformation (rate) is imposed. At low temperatures they behave in a brittle elastic man- 
ner; at high temperatures, much above the glass transition the behavior is more (rubbery) 
(viscoelastic). There is a huge drop in modulus when the temperature is increased above 
the glass transition temperature, indicating a shift in behavior from (glassy) to (rubbery). 
Because of these peculiar mechanical properties of amorphous materials, the linear theory 
of viscoelasticity is unable to model closely the observed response and thus there is a need 
for a non-linear theory of viscoelasticity. 

Unlike crystals, glasses also age, meaning that their state depends on their history. When a 
glass falls out of equilibrium, it evolves over very long time scales. Motivated by the above 
issues and the interesting characteristics of deformation kinetics at very low temperature 
we shall study the deformation kinetics in microscopic domain at low enough temperature 
which is an interesting topic for applications in micro- and nanodomains or the validation in 
using quantum mechanic formulations [13] where the nonlinear constitutive relations should 
be adopted. 

However, real surfaces are rough at the micro- or even at the meso-scale and the role of 
surface roughness has been extensively investigated, and opposite conclusions have been 
reached so far [14-15]. For instance, friction can increase when two opposing surfaces are 
made smoother (this is the case of cold welding of highly polished metals). On the other 
hand, friction increases with roughness when interlocking effects among the asperities come 
into play. This apparent contradiction is due to the effects of length scales, which appear to 
be of crucial importance in this phenomenon. 

Prom the mechanical point of view, a contact problem involves the determination of the 
traction distributions transmitted from one surface to the other, in general involving normal 
pressures and, if friction is present, shear tractions, according to an appropriate set of equal- 
ities and inequalities governing the physics of the contact [16]. When there is friction at the 
contact interface. Coulomb friction behaviour is usually introduced to give the conditions 
necessary to determine the shear traction distribution. Any point in the contact area must 
be either in 'stick', or 'slip' condition, and the tangential tractions must behave accordingly. 
In this paper we shall consider the deformation kinetics of amorphous solid helium at very 
low temperature in micropores which have radius- or transverse-corrugations along the cross- 
section. The glassy matter will be treated as a shear-shinning material. To consider the 
transport of this kind of glass (shear-thinning) material in microdomains, we adopt the 
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verified model initiated by Cagle and Eyring [9] which was used to study the anneahng of 
glass. To obtain the law of annealing of glass for explaining the too rapid anneahng at the 
earliest time, because the relaxation at the beginning was steeper than could be explained 
by the bimolecular law, Cagle and Eyring [9] tried a hyperbolic sine law between the shear 
(strain) rate : F and (large) shear stress : r and obtained the close agreement with experi- 
mental data. This model has sound physical foundation from the thermal activation process 
(Eyring [10] already considered a kind of (quantum) tunneling which relates to the matter 
rearranging by surmounting a potential energy barrier; cf.. Fig. 1) and thus it might also 
resolve the concern raised by Anderson [5] for the thermal noises to the superflow of vortex 
liquid (i.e., the supersolid helium). 

With this model we can associate the (glassy) matter with the momentum transfer between 
neighboring atomic clusters on the microscopic scale and reveals the atomic interaction in 
the relaxation of flow with (viscous) dissipation. 

The outline of this short paper is as follows. Section 2 describes the general physical for- 
mulations of the framework. In this Section, explicit derivations for the glassy deformation 
kinetics are introduced based on a microscopic model proposed by Eyring [10]. The bound- 
ary perturbation technique [17] will be implemented, too. In the third Section, we consider 
the very-low temperature limit of our derived solutions which are highly temperature as 
well as activation energy dependent at rather low temperature regime. Relevant results and 
discussion are given therein. 

II. THEORETICAL FORMULATIONS 

The beginnings of theoretical molecular mechanisms of deformation in amorphous mate- 
rials and glass arc as old as the subject of atomic mechanisms of deformation and yield in 
metals. The first specific molecular mechanism of deformation for amorphous materials and 
glass was published by Eyring [10] and later, Taylor [18] pubhshed the model of an edge 
dislocation to account for the plastic deformation in metals. However, whilst the theory 
of dislocations and crystal defects has become a major stream in the science of sohd state, 
the corresponding effort applied to this problem in amorphous materials must be considered 
rather small by comparison. 

The molecular theory of deformation kinetics came from a different stream of science than 
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that of structure and motion of crystal defects (in particular dislocations) . Its roots stretch 
to the developmental stages of theories of chemical reactions and thermodynamic descrip- 
tion of their temperature dependence, culminating in the key formulation by Arrhenius of 
the equation for reaction rates. By the beginning of this century the concept of activation 
entropy was included in the model, and it was considered that molecules go both in the 
forward direction (product state) and in the backward direction (reactant state). 
The development of statistical mechanics, and later quantum mechanics, led to the concept 
of the potential energy surface. This was a very important step in our modem understanding 
of atomic models of deformation. Eyring's contribution to this subject was the formal devel- 
opment of the transition state theory which provided the basis for deformation kinetics, as 
well as all other thermally activated processes, such as crystallisation, diffusion, polymeri- 
sation, etc. [10] 

The motion of atoms is represented in the configuration space; on the potential surface the 
stable molecules are in the valleys, which are connected by a pass that leads through the 
saddle point. An atom at the saddle point is in the transition (activated) state. Under the 
action of an applied stress the forward velocity of a (plastic) flow unit is the net number 
of times it moves forward, multiplied by the distance it jumps. Eyring proposed a specific 
molecular model of the amorphous structure and a mechanism of deformation kinetics [10]. 
With reference to this idea, this mechanism results in a (shear) strain rate given by 

where 

Vh = A2A3A, Vjn = A2A3A1, 

Ai is the perpendicular distance between two neighboring layers of molecules sliding past 
each other, A is the average distance between equilibrium positions in the direction of motion, 
A2 is the distance between neighboring molecules in this same direction (which may or may 
not equal A), A3 is the molecule to molecule distance in the plane normal to the direction 
of motion, and r is the local applied stress, AE is the activation energy, h is the Planck 
constant, is the Boltzmann constant, T is the temperature, 14 is the activation volume 
for the molecular event [9-10]. The deformation kinetics of the chain is envisaged as the 
propagation of kinks in the molecules into available holes. In order for the motion of the 
kink to result in a plastic flow, it must be raised (energised) into the activated state and 
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pass over the saddle point. This was the earhest molecular theory of yield behaviour in 
amorphous materials, and Eyring presented a theoretical framework which formed the basis 
of many subsequent considerations. 
Solving Eqn. (1) for the force or r, one obtains: 

which in the limit of small {V / B) reduces to Newton's law for viscous deformation kinetics. 
We consider a steady deformation kinetics of the glassy material in a wavy-rough microtube 
of To (in mean-averaged outer radius) with the outer wall being a fixed wavy-rough surface 
: r — To -|- esin(A;^) where e is the amplitude of the (wavy) roughness, and the wave number 
: k — 27r/L (L is the wave length). The schematic is illustrated in Fig. 2. Firstly, this 
material can be expressed as [9-10,14] F = Fq sinh(r/ro), where F is the shear rate, r is the 
shear stress, and 

• _ 2kBTVu -AE 

is a function of temperature with the dimension of the shear rate, 

2kBT 

r. = ^ (4) 

is the referenced (shear) stress, (for small shear stress r ^ tq, the hnear dashpot constitutive 
relation is recovered and to/Fq represents the viscosity of the material). In fact, the force 
balance gives the shear stress at a radius r as r = —{r dp/dz)/2. dp/dz is the pressure 
gradient along the tube-axis or 2;-axis direction. 

Introducing the forcing parameter 11 = —{ro/2TQ)dp/dz then we have F = Fq sinh(nr/ro). 
As the (shear) strain rate is 

du , , 

{u is the rate of deformation (or velocity) in the longitudinal (^-)direction of the microtube), 
after integration, we obtain 

u — Us -\ — z—^ [cosh n — cosh ( — )], (6) 

ii Tq 

here, Ug is the rate of deformation or velocity over the surfaces of the microtube, which is 
determined by the boundary condition. We noticed that a general boundary condition for 
material deformation kinetics over a solid surface was proposed (cf., e.g., [14]) as 

5u^L'Ml-l)-y\ (7) 

■L r. 



where 5u is the rate of deformation (or velocity) jump over the sohd surface, L° is a constant 
shp length and Fc is the critical shear rate at which the shp length diverges. The value of 
Fc is a function of the corrugation of interfacial energy. We remind the readers that this 
expression is based on the assumption of the shear rate over the solid surface being much 
smaller than the critical shear rate of F^. Fg represents the maximum shear rate the material 
can sustain beyond which there is no additional momentum transfer between the wall and 
material-flow molecules. How generic this behavior is and whether there exists a comparable 
scaling for glassy or amorphous materials remain open questions. 

At small pressure gradient, the shear-thinning matter behaves like a Newtonian flow, while 
at high pressure gradient, the shear-thinning matter flows in a plug-flow type. Such a 
behavior is due to the shear thinning of the material, i.e., the higher the shear rate is, the 
smaller is the (plastic) flow resistance [7]. On the microscale, this shear-thinning matter 
can bridge the Newtonian deformation kinetics to that of the pluglike type and offers us a 
mechanistic model to study the deformation kinetics in micro- and even nanodomains using 
the technique of continuum mechanics. 

With the boundary condition from (cf., e.g., [18]), we shall derive the rate of deformation 
(or velocity) field or deformation kinetics along the wavy-rough microtubc below using the 
boundary perturbation technique (cf. [23]) and dimensionless analysis. We firstly select the 
hydrodynamical diameter Lj. to be the characteristic length scale and set 

r'^r/Lr, Ro = ro/Lr, e' = e/L^. (8) 

After this, for simplicity, we drop all the primes. It means, now, r, Ro, Ri, and e become 
dimensionless. The wall is prescribed as r = Ro+e sm{k6), and the presumed fully-developed 
plastic fiow is along the ^-direction (microtube-axis direction). Along the confined (wavy) 
boundaries, we have the strain rate 

^ ^ ^'d/n)^^^ surface' (^) 
where, n means the normal. Let the rate of deformation u be expanded in e : 

u ^ Uq + eui + e^U2 H , (10) 

and on the boundary, we expand u{ro + edr, 9{= 9q)) into 

dr'^ 

u{r, 9)\^ro+edr,eo) = ^^(^0, 9) + e[drUr{ro, 9)] + e^[—Urr{ro, 9)]^ = 



r R TTr 

^coshl^jlfo+^'-^e^^), ro = (11) 

ii /to 

where the subscript means the partial differentiation (say, Ur = du/dr) and 

f 



lis 



I on surface ~ P ) M Ion surface' (-'-2) 
= L°fo[sinhn(l - E^^)-i/2]. (13) 



Now, on the outer wall (cf., e.g., [23]), the (shear) strain rate 
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■ du ^ Vir — Ro — esmCkd)) v/, /im -r 1 

^ = * = ^" ■ |V(r-fi,-.sinU))| = I' + ^V-'(*'')1-'KI(H..*.- 

k 

e— cos{ke)ue\{R„^edr,e)\ = moJro + ^[uiAro + «o.Jbo sin(A;6')- 



k Ik"^ 
-^Uq^\ro cos{ke)] + e^[--— cos^(/c6')mo.|ro +«2.Uo + sin(/c6') + 

7^ Zi f* 



^^io_|R„sin2(^^) - ^cos(A;^)(MiJ,j„ + Mo«Ji?„sin(fc^))] + 0(6^). (14) 

Considering L\ e case, we presume sinh 11 -C Tc/Fo so that we can approximately 

replace [1 - (rosinhn)/re]"^/^ by [1 + Tq sinh n/(2rc)]- With equations (6), (7), (9), (10), 
(11) and (14), using the definition of the (shear) strain rate F, we can derive the rate of 
deformation (or velocity) field up to the second order. The key point is to firstly obtain the 
slip rate of deformation (or velocity) along the wavy boundaries or surfaces. 
After lengthy mathematical manipulations and using (1 — F/Fc)"^^^ pa 1 + r/(2rc), 

= -^^^[cosh(^) - coshn] + Usiip^, (15) 

ii tio 

ui = Fo sm{k6) sinh II + UgHp^, (16) 

we have 

ZL Q if. Zi c if. 

Usiipo + ^^siipi + e^Msap2 + ^(e^) (17) 

where 

lio. = -rosinh(^r), (18) 

JXo 

Uq^, = -ro^ cosh(^r), (19) 



. n 

Mo... = -ro-^sinh(— r), (20) 

k n • n 

= + Mo.. sin(A;6') cos{k6)uog = --^^o cosh(— r) sm{k9), (21) 



and 



Usc^ -^cos'{k9)uo^ + ^uo^^^sm'{ke) = ^to[^cos'{k9)-^sm\ke)]smh{^r). (22) 
Thus, at r = Rq, up to the second order, 

UsHp = u, = L% sinh n(l + ^) + eto sm{ke) [sinh U + cosh H (1 + i^o)] + 

2 rio 2 

sin2(A;^) - cos\ke) + ^ sin2(A;^)] sinh n(l + Ko) + cosh^ n sin2(A;^)}, 

RoLs Ro Ro R-o Tc 

(23) 

where 

Xo = l + (fosinhn)/fe (24) 

From the rate of deformation (or velocity) fields (up to the second order) , we can integrate 
them with respect to the cross-section to get the volume (plastic) flow rate {Q, also up to 
the second order here). 

Q^J^ J u{r,0)rdrd9 ^Qsmooth + ^Qpo+^^Qp2- (25) 



III. RESULTS AND DISCUSSION 

We firstly check the roughness effect upon the shearing characteristics because there are 
no available experimental data and numerical simulations for the same geometric configu- 
ration (microscopic tubes with wavy corrugations in transverse direction). With a scries of 
forcings (due to imposed pressure gradients) : 11 = Ro{—dp/dz)/{2TQ), we can determine 
the enhanced shear rates (dV/dt) due to forcings. Prom equation (5), we have (up to the 
first order) 

^ = ^[sinhn + esin(fc^)^coshn]. (26) 

dt dt Rq 

The calculated results are demonstrated in Figs. 3 and 4. The parameters are fixed below 
(the orientation effect : sin(A;^) is fixed here). To (the mean outer radius) is selected as the 
same as the shp length L° = 100 nm. The amplitude of wavy roughness is e = 0.04, 0.07, 0.1, 
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the Boltzmann constant (/cs) is 1.38 x 10 Joule/°K, and the Planck constant (h) is 
6.626 X 10-34 Joule • s. 

In each panel, the inner curve is the relevant boundary of the tube or the geometric part 
of the presentation. The distance between the inner and corresponding outer curves is the 
calculated physical shear rate : F. We can observe once the temperature (T) changes a little 
from 0.1 ° to 0.15 °, the enhancement of F becomes at least three orders of magnitude (for 
n = 1, the activation energy : 3 x 10""^^ Joules). Even at very low temperature Fig. 4 
gives very large strain rates which are required to to obtain the necessary strain for plastic 
deformation. Thus, the constitutive relations is highly nonlinear at rather low temperature 
regime [9]. It is worth pointing out that the Eyring model requires the interaction between 
atoms in the direction perpendicular to the shearing direction for the momentum transfer. 
This might explain why our result is orientation dependent. The effect of wavy-roughness 
will be significant once the forcing (11) is rather large (the maximum is of the order of mag- 
nitude of e[ntanh(n)/i?o]). 

To be specific, we can illustrate the shear rate (F) with respect to the temperature (T) once 
we calculate Fq as the latter is temperature dependent (but presumed roughness indepen- 
dent here) which could be traced from equation (1). This is shown in Fig. 5. 
Note that, based on the rate-state Eyring model (of stress-biased thermal activation), struc- 
tural rearrangement is associated with a single energy barrier (height) E that is lowered or 
raised linearly by a (shear) yield stress r. If the transition rate is proportional to the plastic 
(shear) strain rate (with a constant ratio : Cq; T — CoRt, Rt is the transition rate in the 
direction aided by stress), we have r = E/V* + (kBT/V*) ln(f /Cqz/o) or 

where V* = V/j is a constant called the activation volume, ks is the Boltzmann constant, T 
is the temperature, i/q is an attempt frequency or transition rate [8,25], and f is the stress 
rate. Normally, the value of V* is associated with a typical volume required for a molecular 
shear rearrangement. Thus, if there is a rather-small (plastic) flow (of the glass) at low 
temperature environment then it could be related to a barrier-overcoming or tunneling for 
shear-thinning matter along the wavy-roughness (geometric valley and peak served as atomic 
potential surfaces) in cylindrical micropores when the wavy-roughness is present. Once the 
geometry-tuned potentials (energy) overcome this barrier, then the tunneling (spontaneous 
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transport) inside wavy-rough cylindrical micropores occurs. 

To examine the behavior of the shear rate at low temperature regime, we calculate Fq and 
r {CqUq — 5 X 10^° s"-*^) with respect to the temperature T and show the results in Fig. 5. 
For a selected activation energy : 5 x 10~^^ Joule or ~ 10^^ eV (a little bit smaller than 
the binding energy of ^He), we can find a sharp decrease of shear rates around T ~ 0.01°K. 
Below this temperature, there might be nearly frictionless transport of glassy matter. Note 
also that, according to Cagle and Eyring [8], V* — 3V6r/2 for certain material during an 
activation event, where V is the deformation volume, ST is the increment of shear strain. 
If we select a (fixed) temperature, say, T = 0.1°K, then from the expression of tq, we can 
obtain the shear stress r corresponding to above forcings (IT) : 

r = Tosinh-^[sinh(n) + esin(A;^)^ cosh(n)]. (28) 

There is no doubt that the orientation effect (9) is also present for deformation kinetics of 
amorphous matter. For illustration (shown in Fig. 6), we only consider the maximum case 
: |sin(fc^)| = 1. The trend of enhancement due to IT (pressure-forcing) and e (roughness) 
is similar to those presented in Figs. 3 and 4. Wc remind the readers that, due to the 
appearance of tq, we fix the temperature to be 0.1°K and the activation volume : 10~^^ m^. 
In fact, as shown in Fig. 6, the calculated (shear) stress (which is directly hnked to the 
resistance of the glassy matter) also shows a sudden decrease around T ~ 0.05°K especially 
for the case of (Coi^o = 2 x 10^ s"^). Here, the activation volume {V* or Vh) is selected as 
0.2 nm^ [19]. Thus, the nearly frictionless transport of the glassy fluid at low temperature 
environment (relevant to the supersolidity, cf. [20]) could be related to a barrier-overcoming 
or tunneling for shear-thinning matter along the wavy-roughness (geometric valley and peak 
served as atomic potential surfaces) in cylindrical micropores when the wavy-roughness is 
present. Once the geometry-tuned potentials (energy) overcome this barrier, then the tun- 
neling (almost frictionless transport) inside wavy-rough cylindrical micropores occurs. 
We also noticed that, as described in [9-10], mechanical loading lowers energy barriers, thus 
facilitating progress over the barrier by random thermal fluctuations. The simplified Eyring 
model approximates the loading dependence of the barrier height as linear. This Eyring 
model, with this linear barrier height dependence on load, has been used over a large frac- 
tion of the last century to describe the response of a wide range of systems and underlies 
modern approaches to sheared glasses. The linear dependence will always correctly describe 
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small changes in the barrier height, since it is simply the first term in the Taylor expansion 
of the barrier height as a function of load. It is thus appropriate when the barrier height 
changes only slightly before the system escapes the local energy minimum. This situation 
occurs at higher temperatures; for example, Newtonian deformation kinetics is obtained in 
the Eyring model in the limit where the system experiences only small changes in the barrier 
height before thermally escaping the energy minimum. As the temperature decreases, larger 
changes in the barrier height occur before the system escapes the energy minimum (giving 
rise to, for example, non-Newtonian deformation kinetics). In this regime, the linear depen- 
dence is not necessarily appropriate, and can lead to inaccurate modelling. This explains 
why wc should adopt the hyperbolic sine law [9-10] to treat the glassy matter. 
To be specific, our results are rather sensitive to the temperature (T) and the activation 
energy. Fig. 7 shows especially the temperature dependence of the forcing parameter (11) 
if dp/dz is prescribed (say, around 6 x 10^° Pa/m) and the activation volume is 0.2 nm^ 
{vo — 100 nm). We can observe that once T increases 11 decreases. F calculated using 
prescribed 11 and using directly T also differs. 

Finally, we present the calculated maximum velocity (unit : m/s) with respect to the tem- 
perature in Fig. 8. Geometric parameters : and the activation volume are the same as 
those in Fig. 7 and the roughness amplitude e = 0.02, O.ObRo- We consider the effect of the 
activation energy : 1.5 x 10^^ and 2 x 10^^ Joule. Around T ~ 0.35°K, the maximum velocity 
(of the glassy matter) either keeps decreasing as the temperature increases for larger activa- 
tion energy or instead increases as the temperature increases for smaller activation energy! 
The results presented in Fig. 8 might be related to the microscopic origin for physical aging 
or effects of thermal history [21] and indeterminate solutions discussed in [22]. The latter 
observation might be related to the argues raised in [23] for the annealing process of solid 
helium at similar low temperature environment if we treat the solid helium to be glassy at 
low temperature regime. On the other hand, we like to remind the readers about the role 
of roughness again. If ^He is isotropic (prepared), and the impurity of ^He distribution was 
added or mixed into ^He quite regularly (equally-distributed). It's possible that the role of 
^He [24] plays under certain selected activation energy and activation volume is qualitatively 
the same as that of the roughness here. In general, as commented in [25] : An understand- 
ing of the NCRI effect may therefore require explicit treatment of the effects of strain on the 
crystal properties • ■ ■ We have used a narrow x-ray beam to study the defects in ^He crystals, 
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and we find crystals contain a microstructure of mosaic regions consistent with small angle 
grain boundaries ■ ■ ■, our presentation here is crucial to the understanding of strain upon the 
sohd hehum at rather low temperature regime once the solid helium is presumed to be glassy 
(cf. [26] for the detailed references) within confined microdomains. To give our approach 
another comparison with the transport data in [3], e.g., with the same temperature 0.4 K as 
in Fig. 4 of [3], once we set the activation energy {A.E) to be 1.0 xlO~^^ Joule (or ~ 10~^ 
eV; cf. Fig. 3 in [27]), the radius to be ~ 1/xm, with other parameters the same as those in 
Fig. 8 above, then we have the maximum velocity being around 100//m/sec which is almost 
the same as that estimated in [3] (where the uniform flow was presumed and the fraction of 
the helium that can flow is 2.5 x 10^^). This is illustrated in Fig. 9. 

IV. CONCLUSIONS 

To conclude in brief, we analytically obtain a class of temperature as well as activation 
energy dependent fields of the rate of deformation for glassy material (like solid hehum) in 
microscopically confined wavy-rough domain at very low temperature regime. The effects 
of wavy corrugation upon the confined deformation kinetics at very-low temperature are 
clearly illustrated. It is found that there exist almost frictionless plastic flow fields for the 
rate of deformation of glassy material inside cylindrical micropores at very low temperature 
once the micropore surface is wavy-rough and the activation energy is prescribed. Once the 
temperature, activation volume, and geometry are fixed, the increase of activation energy 
instead reduces significantly the (maximum) rate of deformation of the glassy matter. The 
critical rate of deformation is proportional to the (referenced) shear rate, the slip length, 
the orientation and the amplitude of the wavy-roughness as illustrated above. Our ap- 
proaches can recover those transport results reported in [3] for the same physical parameters. 
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Fig. 1 Increasing strain causes a local energy minimum to flatten until it disappears 
(removing of energy barrier or quantum-like tunneling) . The structural contribution to 

the shear stress is shear thinning. 
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2. Schematic of a cylindrical micropore, e is the amplitude of small wavy-roughness. 
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Fig. 3. Comparison of the shear rate (F) of glassy matter in smooth microtubes and wavy-rough 
microtubes for A; = 10, e = 0.0,0.04,0.07,0.1 with ro=100 nm. Tc/fo = 10 and L° = Tq. k is the 
wave number and e is the amphtude of the wavy-roughness. T is the temperature. The soHd-hne 

length represents the scale of F = 10s~^. 
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Fig. 4. Comparison of the shear rate (F) of glassy matter in smooth microtubes and wavy-rough 
microtubes for A; = 10, e = 0.0,0.04,0.07,0.1 with ro=100 nm. Tc/Fo = 10 and L° = r^. k is the 
wave number and e is the amphtude of the wavy-roughness. T is the temperature. The solid-hue 
length represents the scale of F = 3000s~^. As the temperature increases a little, F increases 

significantly. 
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Fig. 5 Comparison of calculated shear (strain) rates using an activation energy 5 x 10 
Joule or ~ 10^^ eV. There is a sharp decrease of shear rate around T 0.01°K. 
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Fig. 6 Comparison of calculated (shear) stresses using an activation energy 5 x 10 
Joule or ~ 10"^ eV. There is a sharp decrease of shear stress around T ~ 0.05°K for 
Cq'^o = 2 X 10^ s~^. Below 0.05°K, the transport of glassy matter is nearly frictionless. 
1^0 is an attempt frequency or transition rate [10,19]. 
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Fig. 7 Comparison of calculated shear (strain) rates using an activation energy 10^^^ 
Joule or 10^'^ eV. The temperature (T) dependence of 11 (forcing) is also demonstrated. 
Forcing (11) decreases as the temperature (T) increases. F (mark : diamond) calculated 
using prescribed 11 is different from that (mark : cross) using directly T (temperature). 
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Fig. 8 Comparison of calculated (maximum) velocity (unit : m/s) using two activation 
energies 1.5 x W~^^ and 2 x W^^^ Joule. Around T ~ 0.35°K, the monotonic trend of 
velocity (or deformation rate) bifurcates as the temperature increases. 
ro = 100 nm and e = 0.02, 0.05i2o- 
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Fig. 9 Comparison of calculated (maximum) velocity (unit : m/s) using two activation 
energies 1 x 10~^^ and 1.5 x 10~^^ Joule. Around T ~ 0.25 or 0.35°K, the monotonic trend 

of 

velocity (or deformation rate) bifurcates as the temperature increases. 

To = 1/um and e = 0.02, 0.05i?o. 
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